Poiyeiectroiytes in electric fields: Measuring the dynamical 
effective charge and effective friction 
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We use a coarse-grained molecular dynamics 
model to study the electrophoretic behaviour of 
flexible polyelectrolyte chains. We first charac- 
terize the static properties of the model with 
respect to the chain length, the polyelectrolyte 
concentration, additional salt and the influence 
of an applied external field. Next we inves- 
tigate the dynamic behaviour in the oligomer 
range and compare to data obtained by two 
different experimental methods, namely capil- 
lary electrophoresis and PFG-NMR. We find 
excellent agreement of experiments and sim- 
ulations when hydrodynamic interactions are 
accounted for in the simulations. We then 
present novel estimators for the dynamical ef- 
fective charge during free solution electrophore- 
sis and compare them to static estimators. We 
find complete agreement between the static and 
the dynamic estimators. We further evaluate 
the scaling behaviour of the effective friction of 
the polyelectrolyte-counterion complex with the 
surrounding fluid. We identify a hydrodynamic 
screening length beyond which the friction dur- 
ing electrophoresis is linear depending on the 
chain length resulting in a constant mobility for 
long polyelectrolyte chains. Our results show 
a convincing agreement with experimental data 
and demonstrate that it is possible to model dy- 
namic behaviour of poiyeiectroiytes using coarse 
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grained models, provided they include the ef- 
fects of hydrodynamical interactions. 



1 Introduction 

In order to be able to improve the processes 
involved in current electrophoretic separation 
methods it is a prerequisite to gain a thor- 
ough understanding of the behaviour of poiy- 
eiectroiytes in an externally applied electric 
field. Electrophoretic methods are widely used 
to study poiyeiectroiytes (PEs) such as proteins, 
DNA and syntheticpolymers [l|, 0, Hf. Sev- 
eral theories 0, S, 0| have been developed 
to describe PE electrophoresis and successfully 
described qualitatively the experimentally ob- 
served behaviour of various PEs under bulk con- 
ditions. In a recent publication [8|, we showed 
how coarse-grained molecular dynamics simula- 
tions can extend this theoretical understanding 
on a more microscopic level. In this article, we 
will extend this work by rigorously studying ad- 
ditional effects on the electrophoretic behavior 
such as polyelectrolyte concentration, concen- 
tration of added salt, and variations of the ex- 
ternally applied field. 

We start to characterize the polyelectrolyte's 
static properties on the externally applied field. 
We then continue to investigate the dynamic 
properties, putting emphasis on the hydrody- 
namic interactions between the polyelectrolyte 
and the surrounding solvent. Our results on 
the PE diffusion coefficient and the mobility are 
compared to two experimental data sets and 
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yield very good agreement. Since the mobil- 
ity of the PE can be expressed as the quotient 
of an effective charge and an effective friction 
term, we concentrate on finding a reasonable op- 
erational definition of an effective charge. This 
would then yield a measure of the effective fric- 
tion, and we believe that a direct measurement 
of the effective charge will help to fully un- 
derstand the electrophoretic process. We also 
would like to clarify if there are differences in 
the obtained values of the effective charge, when 
purely dynamical or static estimators are used, 
as has been found for charged colloidal suspen- 
sions [s| . We compare five different ways to esti- 
mate the effective charge of the polyelectrolyte- 
counterion complex, and use the effective charge 
to characterize the effective friction, thereby 
identifying the role of counterions on hydrody- 
namic shielding. 



The paper is organised as follows: in Section [2] 
we present our simulation model and briefly de- 
scribe the different systems investigated. The 
model is characterised in Section [3] by analyzing 
the dependence of static system properties on 
the the system parameters such as chain lengths 
and concentration, effect of added salt, and the 
influence of the external electric field applied 
during electrophoresis. We then determine two 
relevant transport coefficients, the diffusion and 
the electrophoretic mobility, and compare the 
results to recent experimental data (Section Hj). 
To understand the observed results, we intro- 
duce and compare several estimators for the ef- 
fective charge of the polyelectrolyte-counterion- 
complex that is created in solution (Section [5]). 
In Section [6], we use the charge estimate to cal- 
culate the effective friction of the complex dur- 
ing free solution electrophoresis. We identify 
a linear increase of the friction with the chain 
length beyond a screening length. We conclude 
the paper by interpreting the observed results 
and show how this helps to understand the un- 
derlying processes in free solution electrophore- 



2 Model 

We employ molecular dynamics (MD) simula- 
tions using the ESPResSo package [ToJ to study 
the behaviour of linear PEs of different lengths. 
The PEs are modelled by a totally flexible bead- 
spring model in a set of spheres that represent 
the N individual monomers which are connected 
to each other by finitely extensible nonlinear 



elastic (FENE) bonds [11 
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with stiffness k = 30eo, and maximum extension 
R = 1.5cr , and r being the distance between 
the interacting monomers. Additionally, a trun- 
cated Lennard- Jones or WCA potential [l2[ 




r ) + 4 



is used for excluded volume interactions. A cut- 
off value of r c = \/2ao ensures a purely repulsive 
potential. All dissociated counterions and ad- 
ditional salt ions are modelled by appropriately 
charged spheres using the same WCA potential. 
Here, e and a define the energy and length 
scale of the simulations. We use e = k B T, 
i.e. the energy of the system is expressed in 
terms of the thermal energy. The length scale o"o 
defines the size of the monomers. In this model, 
the average bond length between the successive 
monomers in a PE chain is 0.91<To, whereas the 
real distance between two monomers of a fully 
sulfonated polystyrene chain is approximately 
2.5 A. Therefore, loo is mapped to 2.5A/0.91 = 
2.75A in reality. Different polyelectrolytes can 
be mapped by changing cr . Unless mentioned 
otherwise, all observables are expressed in re- 
duced simulation units, and we will not use <Jq 
and eo explicitly from now on. 
All chain monomers carry a negative electric 
charge q = — leo, where eo is the elementary 
charge. For charge neutrality, N monovalent 
counterions of charge +le are added. Where 
mentioned, additional monovalent salt is added 
to the simulation. Full electrostatic interactions 
are calculated with the P3M algorithm using the 



sis. 



implementation of Reference [131 ] . The Bjer 
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rum length l B = ejj/ (47re e r A;B^) = 2.58 in sim- 
ulation units corresponds to 7.1 A, the Bjer- 
rum length in water at room temperature. This 
means that the effect of the surrounding wa- 
ter is modelled implicitly by simply using the 
dielectric properties of water, having a relative 
dielectric constant of e r ~ 80. 
The simulations are carried out under periodic 
boundary conditions in a cubic simulation box. 
The size L of the box is varied to realize a con- 
stant monomer concentration cpe independent 
of chain length as specified. 
We compare two types of MD simulations. In 
the first one, we use simple Langevin equa- 
tions of motions with a velocity dependent dis- 
sipative and a random term in addition to 
the interparticle forces. Together, both addi- 
tional terms implicitly model the effects of a 
solvent surrounding the particles: the dissipa- 
tive force, F D = — T v, provides local fric- 
tion and the non-correlated zero-mean Gaus- 
sian random forces, Fr = —£(£), mimic ther- 
mal kicks (Brownian motion). In order to ful- 
fill the fluctuation-dissipation theorem, dissipa- 
tive and random force have to be coupled to- 
gether: (&(t) (*')> = ^^5^5(1 -t'). This 
approach only offers local particle-fluid interac- 
tions, but destroys long-range hydrodynamic in- 
teractions (HI). 

The second set of simulations include hydrody- 
namics using a Lattice Boltzmann (LB) algo- 
rithm [l4| that is interacting with the MD sim- 
ulations via a frictionally coupling introduced 
by Ahlrichs et al. 15]. The mesoscopic LB 



fluid is described by a velocity field generated 
by discrete momentum distributions on a spa- 
tial grid rather than explicit fluid particles. We 
use an implementation of the D3Q18 model with 
a kinematic viscosity v = 3.0, and a fluid den- 
sity p = 0.864. The resulting fluid has a dy- 
namic viscosity rj = pv = 2.592. The space 
is discretised by a grid with spacing a = 1.0. 
As in the Langevin approach, the particle-fluid 
interaction is realised by a dissipative force de- 
pending on the difference between the particle 
velocity v and the fluid velocity at the particle 
position u: Fr = — rb arc (v — u). Here, the cou- 
pling constant is r bare = 20.0. Additional ran- 
dom fluctuations for particles and fluid act 



thermostat. The interaction between particles 
and fluid conserve total momentum, and this 
algorithm has been shown to yield correct long- 
range hydrodynamic interaction between indi- 



vidual particles [15 



By comparing both types of MD simulations, we 
can characterize the impact of HI on the dynam- 
ics of the system. The individual approaches 
are mapped by setting the Langevin friction pa- 
rameter To = 15.34 to match the single particle 
mobility obtained with the mesoscopic fluid. 
The simulations are carried out with a MD time 
step tmd = 0.01 and LB time step r LB = 0.05. 
After an equilibration time of 10 6 steps, 10 7 
steps are used for generating the data. The 
time-series are analyzed using auto-correlation 
functions to estimate the statistical errors as de- 
tailed in Reference 
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Error bars of the order 
of the symbol size or smaller are omitted in the 
figures. Up to ten independent simulations are 
carried out for each data point, taking between 
one day and three weeks on a single standard 
CPU depending on the chain length N and the 
monomer concentration cpe- 



3 Static properties 

In this section, the proposed coarse-grained 
model is used to quantify the static properties 
of the polyelectrolyte chain and the surround- 
ing counterions. The influence of the individ- 
ual concentrations of the monomers Cpe and the 
counterions cqi as well as the dependence on the 
strength of the external electric field E are in- 
vestigated. 



3.1 Chain scaling 

The polyelectrolyte chain conformations can be 
characterised by the average end-to-end dis- 
tance 



the radius of gyration 
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and by the hydrodynamic radius 



1 

Rh 
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Here, N is the number of chain monomers, the 
position of the z-th monomer, and r cm the cen- 
ter of mass of the polyelectrolyte chain. The an- 
gular brackets (...) indicate an ensemble aver- 
age. All three quantities are expected to exhibit 
a power law scaling -R e , g ,h ~ (N — l) v ', where 
the scaling exponent v depends on the sys- 
tem. For an uncharged polymer with ideal chain 
behaviour v ~ 0.588 (Flory exponent) (l7| . 
whereas for a fully charged polyelectrolyte with- 
out electrostatic screening (no counterions) , we 
expect to find v — 1. 

Figure [TJ shows -R Cjgi h for polyelectrolyte chains 
of different length in the presence of the neu- 
tralizing counterions at a concentration for 
monomers and counterions of Cpe = Cci = 10 
niM and in the absence of an external field. The 
obtained effective scaling exponent v ps 0.85 
for R e and R g for chains up to A^ = 256 in- 
dicates the counterion influence on electrostatic 
screening and on the chain conformations. The 
hydrodynamic radius i?h exhibits a very slow 
asymptotic behaviour which leads to a lower ap- 
parent scaling exponent for short and interme- 
diate chains. In Section 14.11 we will demon- 
strate, however, that this is in perfect agree- 
ment with the measured diffusion coefficients. 
Figure [T] also examplifies that there are no in- 
fluences of the hydrodynamic interactions on all 
static chain properties, as it should be. 
If the polyelectrolyte chain is subject to a strong 
external electric field E, then, depending on the 
magnitude of E, conformational changes can be 
induced 

In Figure [21 we show that the model polyelec- 
trolyte of length N = 40 experiences confor- 
mational changes if the external electric field is 
E, 



crit 



0.2 or larger. According to [18j, Il9l|. E, 



v crit 



depends on the strength of the electrostatic cou- 
pling and on the length N. The approximation 
of Netz E cr i t p 



! Ib/N ps 0.25 yields a value 
comparable to the one found in our simulations. 



The increased end-to-end distance (Figure 2(a) ) 
indicates that the polyelectrolyte chain confor- 
mation is extended. This can also be seen by 



looking at R g in Figure 2(b) Additionally, we 



introduce two new observables 
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and 



R 
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i,o:=y,z 



where 



is the position component in in- 



direction of the i-th particle. For an isotropic 



chain, R 
ure 



R 



i? g yz , which is shown in Fig- 
2(b) for E < 0.1. For electric fields beyond 



g-x 



the threshold, an increase of i? giX and a decrease 
of -Rg, yz is observed. This can be understood 
as an extension of the polyelectrolyte confor- 
mation in x-direction (into the direction of the 
external field) and a compression perpendicu- 
lar to it. This alignment in high electric fields 



has also been studied in Reference 20. There 



it was shown, that for even higher fields and 
stiff polymers, the phase of alignment along the 
electric field is followed by an alignment perpen- 
dicular to it. This transition has also been re- 
cently observed by [2l|. Figure[2]also shows that 



this effect is independent of long range hydro- 
dynamic interactions, which is consistent with 
the initial observations by Netz that were ob- 
tained without the inclusion of hydrodynamic 



18, 19 



interactions 

For the purpose of this paper, we note that for 
electrical fields of E = 0.1 or lower, no confor- 
mational change and no orient at ional order is 
induced. 

Adding salt to a polyelectrolyte solution screens 
electrostatic interactions on a length scale 
known as the Debye length Id- For monovalent 
salt, Id is inversely proportional to the square 
root of the added salt concentration: 



/d 1 = ^4tt/ b (2c s + cci) N A , 

where c s is the concentration of the monova- 
lent salt, cci is the concentration of the poly- 
electrolyte's counterions and is the Avo- 
gadro number. In Figure [31 the screening ef- 
fect of added salt can be seen. We compare 
the scaling of R e and R g for the salt-free case 
(Id ps 16.7) to a solution with c s = 1 M of added 
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salt (Zd ~ 1-2). The additional salt screens the 
electrostatic interactions between the polyelec- 
trolyte monomers and reduces the scaling coef- 
ficient to v — 0.68, which is close to the scal- 
ing for an uncharged polymer with the Flory 
exponent. It remains to be pointed out, that 
the measured scaling exponent is only an effec- 
tive value for the length range investigated. For 
chains with a blob diameter much larger than 
the finite electrostatic screening length Zd, the 
scaling is expected to be equal to the scaling of 
an uncharged polymer, v = 0.588. 
By adding salt, the polyelectrolyte conforma- 
tions become less extended. The effective scal- 
ing exponent v depends on the inverse Debye 
length in the system. We illustrate this in Fig- 
ure IU for no additional salt (small inverse De- 
bye lengths), the observed scaling exponent is 
close to 1 as it is expected for an unscreened 
polyelectrolyte chain. As we add additional salt 
and thereby increase the inverse Debye length, 
the scaling exponent decreases and approaches 
the Flory number. The additional salt screens 
the electrostatic interactions along the polyelec- 
trolyte chain which starts to assume configura- 
tions close to those of an uncharged polymer in 
ideal solvent. 

We note that adding monovalent salt to the so- 
lution, decreases the spatial extension of the 
polyelectrolyte chain and makes the conforma- 
tions more compact. 

3.2 Counterions 

After having investigated the chain conforma- 
tions in some detail, we now take a look at 
the counterion distribution around the poly- 
electrolyte chains. Strongly charged polyelec- 
trolytes attract some of the released counterions 
and effectively reduce their line charge density. 
This phenomenon was described by Manning 
and Oosawa under the term counterion conden- 
This topic has been discussed 



sation 
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sec 



24, 25 and the 



from varying viewpoints 
references therein). 

According to Manning's theory, the distribution 
of counterions around highly charged rodlike 
polyelectrolytes can be described in terms of the 
Manning parameter £ = Zb/&, where Zb = 2.58 



is the Bjerrum length and b = 0.91 is distance 
between charges along the backbone of the poly- 
electrolyte (inverse line charge density), thus 
£ = 2.84. 

Here, we investigate highly charged polyelec- 
trolytes with £ > 1, for which a finite number 
of counterions is always found in close vicinity 
of the polyelectrolyte chain, thus reducing the 
effective charge of the created polyelectrolyte- 
counterion complex. The simple counter ion 
condensation theory predicts the fraction of 
those condensed counterions to be fci = 1 — 1/£, 
but without specifying the actual distance to 
the chain in which those counterions are to be 
found. Following the prediction, the total num- 
ber is 

N CI = (1-1/0 N, (1) 

where N is the length of the polyelectrolyte. 
This leads to a predicted effective charge 

Q eff (iV) = l/^N. 

In Figure we compare this prediction to 
the average number of counterions that we 
find within a distance of 2a to the polyelec- 
trolyte chain. We will present more elaborate 
estimators based on Poisson-Boltzman theory 
in Section [5j For long polyelectrolyte chains 
(N > 100), the measured number of condensed 
counterions approaches the predicted Manning 
value. The distribution of counterions around 
the chain is not influenced by hydrodynamic in- 
teractions, as it should be. 
An external electric field not only couples to 
the polyelectrolyte monomers, but also acts on 
the oppositely charged counterions. Strong elec- 
tric fields are known to reduce the number of 
condensed counterions in the vicinity of the 



chain [18|, 119J. In Figure [6j the number of 
counterions within 2o"o around a polyelectrolyte 
chain with N = 40 monomers is determined for 
different values of the applied electric field. Be- 
yond a critical threshold of E = 0.1, counteri- 
ons are stripped away from the polyelectrolyte 
chain. The onset of this effect coincides with the 
observed extension and alignment of the chain 
(Figure E]). 

Again, we note that for electrical fields of E = 
0.1 or lower, no change to the counterion dis- 
tribution around the polyelectrolyte chain is 
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found. 

4 Transport coefficients 

In this section, we will determine the transport 
coefficients for flexible polyelectrolytes. We will 
briefly review the results from an earlier study 
combining different experimental data sets on 
polystyrene sulfonate (PSS) with simulation re- 
sults [8]. We augment the study by an addi- 
tional investigation about the dynamical and 
static effective charge and friction, respectively, 
to expand the understanding of the microscopic 
processes. 

4.1 Diffusion 

The diffusion coefficient D characterizes the 
thermal motion of the polyelectrolyte. It can 
be obtained from the simulation trajectory of 
the polyelectrolyte chain by recording the center 
of mass mean-square displacement versus time, 
and subsequently measuring the corresponding 
slope according to 



D 



''cm 

(*) - f cm (o)r 

Gt 



(2) 



where r cm is the position of the center of mass, 
and t is the time. The angular brackets (...) in- 
dicate the averaging over many configurations. 
Alternatively, the diffusion coefficient D can be 
obtained from the integration of the velocity 
auto-correlation function of the center of mass 



D 



1 

3 Ju 



(t) -v cm {0))dt. 



(3) 



Here, v cm is the center of mass velocity of the 
polyelectrolyte at a given time. Again, the an- 
gular brackets (...) indicate the averaging over 
many configurations. 

The accuracy of both methods depends on the 
number of statistically independent data sam- 
ples. In Figure we present sample graphs 
to determine the diffusion of a polyelectrolyte 
chain of length N = 32 using Equations [2] and [31 
Since simulations with hydrodynamic interac- 
tions are computationally very demanding, the 
achievable accuracy is limited. The errors are 



determined from the statistical fluctuations of 
the data and the uncertainty in the fit parame- 
ter. 

In Figure 7(a) , the diffusion is obtained from 



a fit to the linear part of the mean square dis- 
placement. This yields a diffusion coefficient of 
D = 0.0045 ± 0.0002 in simulation units. To 
calculate the integral in Equation [31 a fit to 
the slowly decaying long-time tail of the cen- 
ter of mass' velocity auto-correlation function 



has to be obtained as shown in Figure 7(b) 



Here, the theoretical predicted functions are 
used to match the long-time tail: without hy- 
drodynamic interactions, an exponential decay 
of the velocity correlations is expected, whereas 
with hydrodynamic interactions the correlation 
function decays with t~ 3//2 . The figure shows 
the more interesting case with hydrodynamic 
interactions. The combined results of simula- 
tion data and long-time fit are integrated and 
a diffusion coefficient of D = 0.0041 ± 0.0005 is 
obtained, which is in agreement with the corre- 
sponding value obtained from the mean square 
displacement. 

Both methods are strictly equivalent for classi- 
cal systems, but for the remainder of the paper, 
the integral method is used to obtain the diffu- 
sion coefficient as one can use a similar formula- 
tion to obtain the mobility of the polyelectrolyte 
(see Equation [6]), and thus can determine both 
quantities without additional computational ef- 
fort. Specifically, both auto-correlation func- 
tions have to be determined accurately on the 
interval t = [0,100], whereas the root mean 
square displacement has to be determined on 
the interval t = [100, 100000], as shown in Fig- 
ure [7J 

The diffusion of polymers in the presence or ab- 
sence of hydrodynamic interactions is very well 
studied from a theoretical point of view. In gen- 
eral, the following Einstein equation is valid 



k B T 



D 



where T is the so-called friction coefficient of the 
studied object, k B is the Boltzmann factor and 
T is the temperature. 

Without hydrodynamic interactions, one ex- 
pects Rouse behaviour, i.e. a friction coefficient 
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r that linearly depends on the change length N: 



D 



Here, To is the friction coefficient of a single 
monomer of the polymer chain. 
With hydrodynamic interactions, Zimm be- 
haviour is expected. The scaling of D with the 
chain length is no longer proportional to iV" 1 , 
but can be described by the Kirkwood-Zimm 
theory 26|, |27||. Within this theory, the diffu- 



sion coefficient is expected to be 

k h T 



D 



N QTT7]R h ' 



(4) 



where Do is the diffusion coefficient of a single 
monomer of the polymer chain, N is the chain 
length, r\ is the viscosity of the solvent, and 
is the hydrodynamic radius of the polymer. In 
general as pointed out in Section [31 is not 
linear in N, resulting in a scaling different from 
iV -1 for the diffusion coefficient of a polyelec- 
trolyte in the presence of hydrodynamic inter- 
actions. 

In Figure [HJ we compare the normalised dif- 
fusion coefficient for polyelectrolyte chains of 
varying length with and without hydrodynamic 
interactions. When hydrodynamic interactions 
are present, the polyelectrolyte diffusion shows 
a behaviour which can be described by a power 
law scaling D = DoN~ m , where Dq is the 
monomer diffusion coefficient. For our model 
polyelectrolyte, a scaling exponent of m = 
0.63 ± 0.01 is observed. Another way to ob- 
tain the diffusion coefficient from is to insert the 
statically (or dynamically) measured hydrody- 
namic radii and insert them into Equation HI 
Both ways of obtaining the PE diffusion con- 
stant are in perfect agreement, demonstrating 
the applicability of Zimm theory. This is espe- 
cially remarkable, as it shows that the counteri- 
ons, which are not considered in Equation HI do 
not directly influence the diffusion of polyelec- 
trolytes. Instead only the conformations of the 
polyelectrolyte itself, i.e. the internal monomer- 
monomer distances, determine the diffusive be- 
haviour. 

The observed scaling exponent is in good agree- 
ment with experimental results for the diffu- 
sion coefficient of polyelectrolytes as we have 



shown earlier [8|. For the standard polyelec- 
trolyte PSS (polystyrene sulfonate), Bohme and 
Scheler ji^] obtained m = 0.64, whereas Stell- 
wagen et. al. [29[ reported a scaling with m = 
0.617. It is worth mentioning that the latter 
value was obtained in the presence of 50 to 100 
niM additional salt and for PSS chains of upto 
20000 repeat units, therefore a smaller scaling 
exponent can be expected as shown in Section [3) 
When hydrodynamic interactions are switched 
off, the polyelectrolyte chain exhibits Rouse be- 
haviour with a scaling of m — 1.02 ± 0.02, 
demonstrating the importance of hydrodynamic 
interactions for the observed diffusion scaling 
with chain length. 

However, as shown in the previous chapter, the 
static chain properties, including the hydrody- 
namic radius R^ do not depend on the pres- 
ence of hydrodynamic interactions. This allows 
for the correct calculation of the diffusion co- 
efficient even in the absence of hydrodynamic 
interactions by means of Equation HI Having 
said that, the direct measurement of D is only 
possible, when hydrodynamic interactions are 
included in the simulation model. 
For this study, we also lowered the monomer 
concentration of the polyelectrolyte from c m = 
100 mM up to 1 mM. Within this range, no 
change of diffusion coefficients was observed 
(data not shown). 

4.1.1 Influence of the electric field 

When an external electric field is applied to the 
polyelectrolyte solution the diffusive motion in 
the direction of the electric field is mixed with 
the induced directed motion. However, it is 
possible to determine the translational diffusion 
coefficient -Dtrans perpendicular to the electric 
field. Figure M shows that for a reduced elec- 
tric field of E = 0.1 no deviation from the 
diffusive behaviour at vanishing electric field is 
found. This is in-line with the findings of Sec- 
tion [31 i.e., that for small enough electric fields, 
the conformation of the polyelectrolyte chains 
and the surrounding counterions are unchanged, 
such that the measured quantities do not de- 
pend on the applied electric field. 
Experimentally [30], it is also possible to de- 
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termine the diffusion coefficient in the direction 
of the field by separating the diffusive from the 
directed motion. The electric field showed no 
influence on the diffusion coefficient. 



4.2 Electrophoretic mobility 

The second transport coefficient of interest is 
the electrophoretic mobility fi. It characterizes 
the motion of the polyelectrolyte in an external 
electric field. 

In capillary electrophoresis experiments [l|, the 
electrophoretic mobility of the solute is deter- 
mined by 11 = -| = M, where v is the velocity, 
E is the electric field, V is the applied voltage, 
L is the total length of the capillary, I is the mi- 
gration (or effective) length up to the detector, 
and t is the detection time of the solute. 
This method can be directly transferred and ap- 
plied to computer simulations. The external 
electric field E is modelled by a constant force 
acting on the charged particles in the solution. 
This causes a directed motion with a certain 
velocity v. From this one can obtain the elec- 
trophoretic mobility 



long simulations times, in order to be able to ac- 
curately separate the directed motion from the 
thermal fluctuations. 

Alternatively, the electrophoretic mobility can 
be calculated from the following Green-Kubo re- 
lation 

1 f°° 
f M = Wr^Q i (vi(0) ■ v cm (r)) dr, (6) 

where the summation is over all charged parti- 
cles (monomers and counterions) in the system, 
and Vi is their individual velocity and qi their 
charge. Here, v cm is the velocity of the center 
of mass of the polyelectrolyte. This approach 
has been successfully applied in simulations of 



charged colloids [3l|, [32 



While there are several theories 

a, a a a 

that have been used to describe qualitatively 
the experimentally observed electrophoretic be- 
haviour of various PEs, there are still some 
open questions to address. Recent experi- 
ments on strongly charged flexible PEs, such as 
polystyrene sulfonate (PSS) and single-stranded 
DNA (ssDNA) of well defined length have 
shown a characteristic behaviour for the short 



v 

E' 



(5) 



Netz [18|, [19J studied the behaviour of flexi- 
ble polyelectrolytes in strong electric fields in 
the absence of hydrodynamic interactions. He 
showed, that the electrophoretic mobility de- 
termined by Equation for large electric fields 
depends on the magnitude of the electric field. 
This effect is attributed to the polarization and 
the following removal of the counterion cloud 
surrounding the polyelectrolyte. Below a crit- 
ical value, the mobility is not affected by the 
electric field, and the system is in the linear re- 
sponse regime. 

Experimentally used electric fields are usually 
below 1 kV/cm, which corresponds to a reduced 
field strength of E = 0.001, and as such are far 
below this critical threshold. Thus we will fo- 
cus on the behaviour of polyelectrolytes in weak 
electric fields. In Section [3j we identified the 
threshold to be E CTit « 0.2. For simulations, the 
usage of a weak electrical driving force requires by Frank and Winkler [21 



chain free-solution mobility /i [29|, |33|, |3J, [35 
after an initial increase of the mobility with in- 
creasing length, fi passes through a maximum, 
and then decreases towards a constant mobility 
for long chains. 

The increase of \x for short chains and the 
long-chain limit constant \x can be explained 
within the theoretical approaches, but the ori- 
gins for the maximum for intermediate chains 
have not been accounted for yet. In [§], we 
showed that the experimentally observed be- 
haviour can be simulated using a coarse-grained 
MD model. We also showed that the maximum 
is due to the complex hydrodynamic interac- 
tions between the polyelectrolyte and its coun- 
terions and the solute. 

In Figure [TQl we illustrate that the maximum 
in the electrophoretic mobility can only be re- 
produced when hydrodynamic interactions are 
properly accounted for. The neglect of hydro- 
dynamic interactions leads to a decreasing elec- 
trophoretic mobility for short chains. This ob- 
servation was also made in a recent publication 
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In addition to the measured mobilities, Fig- 
ure [10] includes a prediction for the mobility 
without hydrodynamic interactions based on a 
local force balance: without hydrodynamic in- 
teractions, every particle of the polyelectrolyte- 
counterion complex is subject to the same fric- 
tional force insolvent = ~^o v that counterbal- 
ances the electric driving force Fpieid = q%E. 
From this we can easily obtain the following 
expression for the electrophoretic mobility in 
absence of hydrodynamic interactions, where 
fiQ = 1 /To is the mobility of a single monomer, 
N is the length of the polyelectrolyte, and iVci is 
the number of condensed counterions that move 
with the polyelectrolyte: 



ii_ = N-Ng 

Ho N + N C1 



(7) 



In order to utilize Equation [7] we need the num- 
ber of condensed counterions. These we obtain 
Nqi by counting the average number of counter- 
ions found within 2<7o of the chain. Additionally, 
we can substitute Nqi from Eq. [Hand obtain the 
Manning prediction for the mobility 



JL 



1 



2£ 



0.2. 



(8) 



Figure [JO] shows excellent agreement of of Equa- 
tion [7] and also demonstrates, that the Manning 
mobility is approached nicely for long chains. 
We conclude that the local force picture suc- 
cessfully describes the observed behaviour in ab- 
sence of hydrodynamic interactions, but quali- 
tatively fails to describe the mobility for shor 
chains in any real experiment, where hydrody- 
namic interactions are obviously present. 

4.2.1 Influence of the electric field 

The effect of high electric fields on the mo- 
bility of polyelectrolytes has been investigated 
theoretically before. Below the critical field 
strength, i.e. in the linear response regime, the 
electrophoretic mobility /i is independent of the 
applied electric field. In Figure [TTJ, we com- 
pare /i for a polyelectrolyte chain at c m ob- 
tained at zero field via Equation [6] to the mo- 
bility at finite fields, E = 0.05 and 0.1, via 
Equation The measured values agree within 



their displayed precision. For comparison, all 
data sets were obtained using the same com- 
putational effort (same number of simulation 
steps). For short chains, the Green-Kubo based 
method yields more accurate results, whereas 
longer chains can be simulated at equal accu- 
racy by an external applied field. For short 
chains and weak fields, the thermal (Brownian) 
motion dominates the directed electrophoretic 
motion, which decreases the accuracy of the ob- 
tained values. 

Based on this observation, we conclude that 
the Green-Kubo method to measure the elec- 
trophoretic mobility of polyelectrolytes in solu- 
tion has several advantages. Firstly and most 
importantly, it guarantees the measurement of 
free-solution mobilities in the linear regime, 
which are comparable to experimental measure- 
ments at standard field strengths. The alterna- 
tive method uses fields that are about a factor 
100 higher than the experimentally used ones 
and are close to the critical value at which static 
and dynamic properties of the polyelectrolyte 
are significantly changed. Secondly, the com- 
putational effort needed to achieve a given ac- 
curacy for short polyelectrolytes chains (N < 
10) is up to 50% smaller than with the direct 
method. And last but not least, the trajectories 
simulated at zero field can be used to determine 
the electrophoretic mobility and the diffusion 
(using Equation [3]) at the same time, addition- 
ally reducing the computational cost of such a 
study. 

4.2.2 Influence of the monomer concen- 
tration 

It is well known, that the free- solution mobility 
of polyelectrolytes is depending on the salt con- 



centration of the solution [34(. With increasing 
additional salt, the maximal free-solution mo- 
bility decreases due to increased counterion con- 
densation. We point out that the long chain mo- 
bility also shows a dependence on the monomer 
concentration in the absence of additional salt. 
This effect can be mainly attributed to the 
electrostatic screening mediated bu the chain 
monomers and counterions (self-screening). 
Figure [J2] shows that not only the long chain 
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limit is influenced by increased electrostatic 
screening, but also the short chain behaviour. 
The maximum is significantly reduced for higher 
monomer concentrations. From a value higher 
or equal to c m = 100 mM it is completely re- 
duced, and the behaviour in the presence of hy- 
drodynamic interactions is similar to one with- 
out hydrodynamic interactions as seen in Fig- 
ure [lOl The increased electrostatic screening 
caused by a higher concentration (resulting in 
a shorter Debye length) suppresses the short- 
range hydrodynamic interactions that are essen- 
tial for the formation of the maximum. 
This is reflected in Figure [T2J where we show 
the degree of polymerization with the maximum 
mobility iV max versus the Debye length Ad- The 
lower the Debye length, i.e. the higher the elec- 
trostatic screening, the more the maximum is 
shifted to shorter chains. This correlation has 
not been investigated previously. In the remain- 
der of this paper, we will analyze the importance 
of the interplay of the electrostatic screening 
with the hydrodynamic interactions. 



5 Charge estimators 

In this section, we present five different esti- 
mators to measure the effective charge Q c s of 
the polyelectrolyte counterion complex. The 
first two are based on static ion distributions, 
whereas the latter three are based on dynamic 
quantities that can be computed with or with- 
out hydrodynamic interactions included. We 
will compare the practicability and accuracy of 
the estimators and discuss the obtained results 
in the context of polyelectrolyte electrophore- 
sis. At the same time we can see if any differ- 
ence can be observed between a static effective 
charge and a dynamic effective charge @, l3l|. 

5.1 Primitive estimate 

A simple method of estimating the effective 
charge has been used in Figure for describ- 
ing the static properties of the polyelectrolyte 
counterion complex: 



where Na(d < do) is the average number of 
counterions that can be found within a distance 
d around the polyelectrolyte, where d is distance 
to the closest monomer. The threshold do is 
usually chosen to be d = 2a . 
In 



Section I3.2[ we showed that the effective 
charge estimated in this way is not affected by 
hydrodynamic interactions. 
This method is computationally inexpensive 
and as seen in Figure \5\ yields a reasonable esti- 
mate for the effective charge. The draw back is 
that the threshold do is arbitrarily defined which 
limits the usefulness of this estimator. 

5.2 Inflection criterion 

A more advanced method uses the inflection cri- 
terion to estimate the threshold of counterion 



condensation [36|, |37|, |38( that has its basis in 
Poisson-Boltzmann theory. It has been shown 
that the position d c of the inflection point (zero 
of the second derivative) for the integrated ion 
distribution with respect to the logarithmic dis- 
tance from the closest chain monomer yields a 
cutoff that accurately separates free from bound 
counterions in the case of a rod-like polyelec- 
trolyte in cylindrical geometry. The applicabil- 
ity to polyelectrolytes will be shown here. 
The integrated ion distribution is defined as 



1(d) 



p cl (r)dr = /ci(r < d), 



(10) 



= JV PE - N CI {d < d ) 



(9) 



where pa is the normalised density of counteri- 
ons at a distance r to the closest chain monomer, 
and fci(r < d) is the fraction of counterions 
found at a distance closer than d. The total 
number of condensed counterions up to this dis- 
tance given by multiplying this number with 
Nci- 

In order to apply the inflection point criterion 
to the polyelectrolyte model, the integrated ion 
distribution has to be measured and plotted log- 
arithmically. Figure dH shows the result for a 
chain of N = 32 monomers. Additionally, the 
first and second derivative obtained numerically 
are plotted. The cutoff value d c indicated by the 
inflection point and the associated fraction of 
condensed counterions fa can be directly read 
off Figure QU 
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The effective charge of the polyelectrolyte using 
the inflection criterion is then given by 



=N PE -N ci {d<d c ). 



The cutoff d c is not a fixed parameter anymore, 
but has to be determined for every chain length. 

5.3 Langevin model 



In Section 14.21 we used the local force picture 
and derived Equation [7] for the electrophoretic 
mobility in absence of hydrodynamic interac- 
tions. Similarly, we can derive an expression 
for the effective charge of the polyelectrolyte in 
absence of hydrodynamic interactions based on 
the measured mobility fi of the PE-counterion 
complex. 



o (3) 

Veff 



N 



PE 



1 



l + /xlY 



(12) 



Equation [12] offers a way to determine the ef- 
fective charge of the polyelectrolyte counterion 
complex based purely on the measured mobil- 
ity in the absence of hydrodynamic interactions 
without the need of free parameter and at rea- 
sonable computational costs. The accuracy of 
this estimator is limited by the accuracy of the 
measured electrophoretic mobility /i, and it is 
strictly valid only for the Langevin model with- 
out HI. 

5.4 Ion diffusion 

The local force picture used for the Langevin 
model assumes that a finite amount of counte- 
rions is bound to the polyelectrolyte while the 
remaining ions can move freely. We will now 
motivate another charge estimator that is espe- 
cially promising as it can also be directly ap- 
plied in experimental setups, since all necessary 
quantities can be measured experimentally. 
The iVci bound counterions are expected to dif- 
fuse together with the polyelectrolyte with a 
diffusion coefficient D PE as determined in Sec- 
tion H]TJ Likewise, the N — Nqi free counterions 
will diffuse with a different diffusion coefficient 
-Dq- If one measures the ion diffusion coefficient 



Dqi for all N ions in such a system, the mea- 
sured quantity will be the weighted average of 
the bound and the free diffusion coefficient: 



D 



ci 



iVciPpE + (N- N cl )D 
N 



From this we derive a novel estimator for the 
effective charge 



o (4) 

Vcff 



Nil 



Dn — D 



ci 



Dn — D 



PE 



(13) 



that has the advantage of only including quan- 
tities that are experimentally accessible and as 
such can be used to directly measure the effec- 
tive charge in experiments. 
In Figure [IB] we compare the results of esti- 
mator QefF with and without hydrodynamic in- 
teractions and that within the accuracy of the 
method, no difference can be observed, even 
though the individual diffusion coefficients are 
different for both types of simulations. Equa- 
tion [13] is valid for simulations independent of 
the presence or absence of hydrodynamic inter- 
actions and has the computational complexity 
of determining the diffusion coefficients (see Sec- 
tion OJ). 

This estimator is especially interesting, as it 
only involves quantities that are directly obtain- 
able from experiments. As such, it might pro- 
vide a great opportunity to measure the effec- 
tive charge of macromolecules. To our knowl- 
edge, this method has not been suggested be- 
fore. The derivation of equation [T3l assumes the 
absence of additional salt, i.e. c s = 0, but can 
be modified to account for a finite amount of 
additional salt. A similar analysis can be done 
using the free and the bound mobility, but these 
quantities are experimentally more difficult to 
obtain. 

5.5 Co- moving counterions 

The last estimator for the effective charge we in- 
troduce in this article is based on directly deter- 
mining the counterions that are co-moving with 
the polyelectrolyte during free-solution elec- 
trophoresis similar to the method used in [39, 
3. 
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When applying an external electric field E, the 
polyelectrolyte moves with a velocity t>PE = 
[j>peE in the direction of the electric field. Co- 
moving counterions move with the same velocity 
in the same direction, whereas free counterions 
move with a velocity t>ci,o — f^oE into the op- 
posite direction. Here, no is the mobility of free 
counterions. 

In Figure [TU we plot the velocity v ci of the 
ions versus the distance d to the center of mass 
of the polyelectrolyte. It can be seen, that ions 
close to the center of mass are co-moving, i.e. 
^ci — ^PE) and ions far away from the chain are 
indeed freely moving with a velocity vci = vqi,o- 
For an intermediate value of d , the vqi = 0, i.e. 
the ions are not moving relative to the electric 
due to the interaction with the polyelectrolyte. 
This distance do separates the co-moving coun- 
terions from the free-moving ones. Similar to 
Qeff (Eq. ED, the effective charge is obtained 
by averaging the number of counterions found 
within a distance do of the center of mass of the 
chain. 



o (5) 

Veflf 



iVpE - N C i(d < d ) 



(14) 



Note, that for this estimator d is the distance of 
the counterion to the center of mass of the poly- 
electrolyte, not to the closest polyelectrolyte 
monomer. 

The results shown in Figure [H] indicate a small 
influence of the hydrodynamic interactions on 
the value of do, but it is important to note that 
the ion density in this region is very low result- 
ing in almost identical values of the in both 
cases. 

The threshold used for this estimator is not pre- 
defined, but has to be determined from the sim- 
ulation and generally is a function of the chain 
length. Since vqi has to be determined for each 
distance d to the center of mass of the polyelec- 
trolyte requiring high statistics. 

5.6 Comparison 

Figure [17] compares the results of all presented 



charge estimators to Qgff ■ For short chains 
(N < 4), the estimators agree and coincide 
with the bare, unscreened charge of the poly- 



j(5) 



electrolyte. In this regime, no counterion con- 
densation is observed. 

For intermediate and long chains (N > 4), the 
effective charge is reduced as it deviates from 
the bare charge and tends towards the Man- 
ning prediction. However, it probably will never 
reach it completely as the polyelectrolyte starts 
to assume a coiled conformation that is not ac- 
counted for in condensation theory. Here, the 
condensation parameter for the polyelectrolyte 
system is £ = 2.84. In this regime, the simple es- 
timator Qeff measures a higher effective charge, 
i.e., not all condensed counterions that are cor- 
rectly included in the other estimators are taken 
into account. This effect is not observed if the 
cutoff value is chosen using the inflection crite- 
rion. 

All estimators show no or little influence of 
the hydrodynamic interactions on the effec- 
tive charge of the polyelectrolyte counterion 
complex. This independence on the hydrody- 
namic interactions has been recently observed 
for highly charged colloid using a similar simu- 



lation approach 39 . 

Furthermore, no difference between the static 
estimate, Q c g and the dynamic estimates Q eS 
for the effect charge are observed. This is 
especially remarkable as it is an open ques- 
tion if there is a difference between the static 
and the dynamic effective charge. For the 
case of charged colloids this seems to be the 
case P, |3l[, but the results of this work show 
that for strongly charged linear polyelectrolytes 
both quantities are identical and one does not 
have to differentiate between a statically and 
dynamically renormalized charge. 
Two of these estimators are particularly promis- 
ing. The charge estimator based on the 

(3) 

Langevin model, Q eS , has the advantage that 
it is only necessary to determine one dynamic 
quantity, namely the electrophoretic mobility. 
All other parameters are directly given by the 
model. In this way one can obtain an accu- 
rate estimate without much computational ef- 
fort. The ion-diffusion estimator Q^g, on the 
other hand, requires the calculation of three dy- 
namic quantities, which is slightly more costly 
in simulation but has the great advantage of be- 
ing directly transferable to experiments, since 
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all needed quantities, i.e., the free ion diffusion, 
the bound ion diffusion, and the chain diffusion 
can be measured by standard techniques. To 
the author's knowledge, this is the first time this 
method is proposed to determine the effective 
charge of polyelectrolytes during electrophore- 
sis. 

6 Effective friction 

We will now use the effective charge obtained 
in the previous section to quantify the effective 
friction of the polyelectrolyte-counterion com- 
plex. As shown in Figure [TTJ, the estimators 
QS to Qes Yield the same effective charge, we 

(3) 

thus only employ Q e g to determine the effective 
friction in the following equation: 



Q 



elf 



eff 



(15) 



Here, we use \i as obtained in Section 14.21 
From the result displayed in Figure [HI we can 
see that for short chains the effective friction de- 
termined by Equation [15] is in agreement with 
the hydrodynamic polymer friction T = 1/D 
that follows the Einstein relation. For longer 
chains, on the other hand, a significant devia- 
tion is observed. The effective friction no longer 
follows the N~ m behaviour of the Einstein rela- 
tion (with m pa 0.63), but instead tends towards 
a linear increase in N. The length scale separat- 
ing both regimes is of the order of the Debye 
length in the system. 

From this we infer that the increased effective 
friction is caused by the counterions close to 
the polyelectrolyte chain, which destroy long- 
range hydrodynamic interactions between dis- 
tant parts of the polyelectrolyte. In this way 
they also increase the hydrodynamic screening, 
and effectively decouple the chain monomers 
from each other. An similar conclusion is 
obtained in an independent study by Fischer 
et. al. 
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The effect is strongly related to the 
counterion density in the vicinity of the chain 
and thus depends on the Debye length in the 
system. 

It becomes even more evident when looking at 
the effective friction per monomer, Figure [TJB 



Initially, for short chains, the effective friction 
decreases strongly with increasing chain length. 
The monomers move together through the sol- 
vent and can shield each other from the flow and 
in this way reducing the effective friction with 
the solvent. For longer chains, this reduction of 
the effective friction becomes less effective be- 
cause the counterions that associate with the 
polyelectrolyte influence the solvent flow around 
the polyelectrolyte effectively canceling the ben- 
eficial shielding effects. On a length scale that 
is comparable to the Debye length in the sys- 
tem, different parts of the polyelectrolyte be- 
come hydrodynamically decoupled and the ef- 
fective friction per monomer does not depend 
on the length of the polyelectrolyte any more. 
The effective charge per monomer shows a dif- 
ferent behaviour. Short chains do not have 
any bound counterions and show their bare 
charge of le per monomer. With increasing 
length, they can attract counterions which are 
then co-moving with the polyelectrolyte and 
reducing its effective charge. Again the con- 
centration of counterions around the polyelec- 
trolyte plays an important role and the rele- 
vant length-scale can be compared to the De- 
bye length. The similarities in length scales of 
the Debye length to the hydrodynamic screening 
length have been treated analytically previously 



in Refs. 42, 43, 44 



The combined behaviour of effective friction and 
effective charge leads to the observed length de- 
pendent mobility for short chains and the con- 
stant mobility for long chains is illustrated in 
Figure [20j Initially, very few counterions are 
attracted to the polyelectrolyte and the rise in 
mobility is due to the hydrodynamic interac- 
tions between the chain monomers, which are 
in an extended conformation. Then, at inter- 
mediate chain lengths, the counterion conden- 
sation increases, reducing the effective charge 
of the polyelectrolyte, which leads to the ob- 
served maximum. For longer chains, the counte- 
rions shield the electrostatic and hydrodynamic 
interactions between the chain monomers and 
thereby induce a transistion to less extended 
conformations, that for even longer chains (not 
shown here) become globular and coil-like. The 
screening of hydrodynamic interactions together 
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with the conformational change increases the ef- 
fective friction of the polyelectrolyte, leading to 
the observed decrease in the electrophoretic mo- 
bility and the long-chain plateau value. 
More specifically, the mobility maximum, which 
is observed for flexible polyelectrolytes, is due 
to the efficient shielding between monomers for 
short chains which reduces the effective friction. 
Stiffer polyelectrolytes experience a stronger 
friction since the monomers in rod-like con- 
formations can not shield each other that ef- 
ficiently. This is the reason, why the maxi- 
mum in the mobility for intermediate chains 
is only observable for flexible or semi-flexible 
polyelectrolytes such as PSS or single-stranded 
DNA but not for the more rigid double-stranded 
DNA. On the other hand, the decrease of the 
effective charge by counterion attraction is de- 
pending on the linear charge density of the 
polyelectrolyte. The higher the polyelectrolyte 
is charge, the more the effective charge is re- 
duced by co-moving counterions. Additionally, 
we showed that the counterions also increase the 
effective friction. Both effects work in the same 
direction and cause the maximum to be shifted 
to shorter chains (in case of PSS) or disappear 
completely (for double-stranded DNA). 



7 Conclusion 

We studied in detail the electrophoretic be- 
haviour of flexible polyelectrolyte chains by 
means of a coarse-grained molecular dynamics 
model. The static chain properties have been 
investigated for various values of the external 
electric field, and for different salt concentra- 
tions. The static properties do not depend on 
the presence or absence of hydrodynamic in- 
teractions. Our results show the expected be- 
haviour of polyelectrolyte chains: the scaling 
of R e , g ,h for short chains exhibits an effective 
exponent that lies between the infinite dilution 
exponent, v = 1, and the Flory exponent for 
uncharged chains, v = 0.588, and v decreases 
with increasing salt concentration due to elec- 
trostatic screening. 

We tested the influence of the strength of an 
applied external electric field on the chain con- 



formations and the counterion cloud. Below a 
critical value £7 C rit, no dependence on the field 
strength is observed, which is in agreement with 
linear response theory. An applied electric field 
larger than a critical value of E CT & ~ 0.2 leads 
to extended chain conformations and to a sig- 
nificant loss of condensed counterions. Our re- 
sults agree with the theoretical predictions by 
R. R. Netz 0, Gjj, which, 

0.25. 



for our system, 

yielded a value of E c 
When simulating the dynamic transport prop- 
erties, such as diffusion coefficient and elec- 
trophoretic mobility, long-range hydrodynamic 
interactions become crucial. Only when they 
are included, the experimentally observed be- 
haviour is reproduced and excellent agreement 
between experiments and simulations is found 
for our PSS model system. Our results demon- 
strate convincingly that it is possible to model 
quant it avively the dynamic behaviour of poly- 
electrolytes using coarse grained models. 
To investigate the dynamics even further, 
we present five different approaches to esti- 
mate the dynamical effective charge of the 
polyelectrolyte-counterion complex during free- 
solution electrophoresis. The charges calculated 
by the estimators are not influenced by hydro- 
dynamic interactions, nor by the fact if they 
were obtained from static or dynamic observ- 
ables. Thus we find no difference between a 
static and a dynamic effective charge. In ad- 
dition, it is possible to use the estimators in 
simulations without hydrodynamic interactions 
that are computationally inexpensive. 
For the charge estimate based on the Langevin 

(2) 

model, Qeff , it is only necessary to determine 
one dynamical quantity, the electrophoretic mo- 
bility. All other parameter are directly given by 
the model. In this way one can obtain an accu- 
rate estimate without too much computational 
effort. The ion-diffusion estimator on the 

other hand, requires the calculation of three dy- 
namical quantities, which is slightly more costly 
in simulation but has the great advantage of be- 
ing directly transferable to experiments since all 
needed quantities, i.e. the free ion diffusion, the 
bound ion diffusion and the chain diffusion can 
be measured by standard techniques. To our 
knowledge, this is the first time, this method 
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is proposed to determine the effective charge of 
polyelectrolytes during electrophoresis. 
Using theses estimators, we determined the 
length dependence of the effective charge, and 
combined it with the measurements of the elec- 
trophoretic mobility to obtain the effective fric- 
tion of the polyelectrolyte. The results indicate 
that this effective friction during electrophore- 
sis is different from the hydrodynamic friction 
for a single polyelectrolyte chain obtained from 
diffusion measurements. We attribute this dif- 
ference to the contribution by the co-moving 
counterions, which cause a shielding of the hy- 
drodynamic interactions. We identified a hy- 
drodynamic screening length beyond which the 
effective friction approaches a constant value 
per monomer which - together with the con- 
stant value per monomer of the effective charge 
- leads to the well-known and observed con- 
stant electrophoretic mobility for long polyelec- 
trolyte chains. We showed that this hydrody- 
namic screening range is comparable to the De- 
bye length for electrostatic screening with the 
system. 

The results of this study provide an in-depth 
understanding of the microscopic processes that 
govern the macroscopic behaviour of charged 
polyelectrolytes in free solution. The charac- 
terisation of the dynamical effective friction of 
macromolecules based on the direct measure- 
ment of the effective charge seems to be a 
promising tool to investigate and understand 
related systems, such as the dynamics of poly- 
electrolytes slowed down by additional hydrody- 
namic drag tags (ELFSE). A study of this topic 
is currently under preparation. 
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Fig. 1 Both the end-to-end distance R e and 
the radius of gyration R g at vanish- 
ing external field for a monomer con- 
centration of cpe = cci = 10 mM 
exhibit an effective scaling exponent 
v = 0.85. The hydrodynamic radius 
i?h shows a different behaviour for 
short chains and only slowly reaches 
the asymptotic scaling. The static 
chain properties are not influenced 
by hydrodynamic interactions (HI) 
and fully agree with the ones ob- 
tained without HI. 

Fig. 2 (a) The end-to-end distance R e and 
(b) the radius of gyration radius R g 
for a chain of A = 40 monomers 
differ from the value at zero exter- 
nal field E (dashed line), if a crit- 
ical value is reached. Beyond this 
threshold, the polyelectrolyte chain 
is extended (increased R e and R s ) 
and starts to align with the field di- 
rected in x-direction. This effect is 
independent of hydrodynamic inter- 
actions (HI). 

Fig. 3 The end-to-end distance R e and the 
radius of gyration R g at zero exter- 
nal field and for a monomer con- 
centration cpe = 10 mM and no 
added salt (c s = M) exhibit an 
effective scaling exponent v = 0.85 
(dashed line). With c s = 1 M mono- 
valent salt added, the scaling expo- 
nent drops down to v — 0.68 (dotted 
line) . 

Fig. 4 The effective scaling factor v shows 
strong dependence of the Debye 
length. With out added salt 
(i.e. large Debye lengths) v is close 
to 1, whereas for high salt concen- 
trations (i.e. small Debye lengths) 
approaches the Flory number, the 
value of an uncharged polymer. 

Fig. 5 As a simple estimate of the number 
of condensed counterions A CI all ions 
within a distance of 2<7o around the 



PE chain are summed up. For long 
polyelectrolyte chains, the value pre- 
dicted by Equation [1] is approached. 
The static counterion distributions 
are not influenced by hydrodynamic 
interactions (HI). 

Fig. 6 The number of condensed counteri- 
ons iVci for a polyelectrolyte chain 
of A = 40 depends on the applied 
external electrical field. Above a 
threshold of E = 0.1, counterions 
start to be stripped from the PE 
chain and A*ci differs from the value 
at zero field (dashed line). This ef- 
fect is independent of hydrodynamic 
interactions (HI). 

Fig. 7 The diffusion coefficient D of a poly- 
electrolyte chain of length A = 32 is 
determined, (a) A fit (dashed line) 
to the linear part of the mean square 
displacement yields a diffusion coef- 
ficient of D — 0.0045 ± 0.0002 via 
Equation [3 (b) By fitting a r 3/2 
power law to the long-time tail of the 
center of mass auto-correlation func- 
tion and using Equation [3] we obtain 
D = 0.0041 ±0.0005. 

Fig. 8 The diffusion of a polyelectrolyte 
chain of length A at a monomer con- 
centration c m = lOmM, normalised 
by the monomer diffusion D , shows 
the influence of hydrodynamic inter- 
actions (HI). With HI a scaling ex- 
ponent of m = 0.63±0.01 is obtained 
(solid line), whereas without HI, the 
slope is m = 1.02 ± 0.02 (dashed 
line) . The diffusion in presence of HI 
agrees with the values obtained from 
Equation HI where R^ is determined 
from the simulation (triangles). 

Fig. 9 The translational diffusion coeffi- 
cient measured at a reduced electric 
field E = 0.1 (circles) is in agree- 
ment with the measurements at zero 
electric field (solid line) 
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Fig. 10 



Fig. 16 



Fig. 11 



Fig. 12 



Fig. 13 



Fig. 14 



Fig. 15 



Fig. 17 



The normalised electrophoretic mo- 
bility at a monomer concentration 
of c m = 5 mM with hydrodynamic 
interactions (HI) reproduces the 
experimentally observed behaviour 
and shows a maximum for intermedi- 
ate chains. Without HI the measures 
mobility strongly deviates from this 
showing a decrease with increasing 
chain length N. This behaviour can 
be explained by Equation [7] in the 
local force picture (dashed line) and 
approaches the limiting value of the 
Manning prediction (Eq. [8]). 

Electrophoretic mobility \x of poly- 
electrolyte chains at monomer con- 
centration c m = 5 mM measured 
without electric field (circles) and 
with field at E = 0.05 and 0.1 (tri- 
angles) with the corresponding error 
bar. (Data sets have been shifted to 
increase readability.) 

The normalised electrophoretic mo- 
bility /i//io shows a strong depen- 
dence on the monomer concentration 
c m . For dilute systems, the experi- 
mental behaviour is recovered, while 
at high monomer concentrations, the 
shape is significantly altered. 

The approximate position of the 
maximum iV max changes with the 
Debye length Ad- 

The integrated ion distribution 1(d) 
for a chain of iV = 32 monomers 
shows a clear inflection point identi- 
fied by the zero of the second deriva- 
tive with respect to \n(d). The value 
of 1(d) at the inflection point speci- 
fies the fraction of condensed coun- 
ter ions fci- 



Estimated charge Qeff °f polyelec- 
trolyte chains of varying length N 
using Equation [13] with and without Fig. 20 
hydrodynamic interactions (HI). 



Fig. 18 



Fig. 19 



The ion velocity distribution vci(d) 
for a polyelectrolyte chain of length 
N = 32 shows the transition from 
co-moving counterions with a veloc- 
ity similar to vpe (dashed line and 
dotted line indicate vpe with and 
without hydrodynamic interactions 
(HI)) to the free ion velocity vci,o 
(dash-dot line). Vci(d = d ) = de- 
fines the threshold between co- and 
counter moving ions and is used to 
determine the effective charge of the 
polyelectrolyte counterion complex 
in Equation [TJ~ 



The effective charge of polyelec- 
trolyte chains of length N using 
the estimators to . The 
dotted line indicates the bare, un- 
screened charge of the polyelec- 
trolyte, whereas the solid line shows 
a prediction based on counterion 
condensation theory, with £ = 2.84 
being the condensation parameter. 
Whereas estimators to 
agree over the range of lengths N 
tested, Qeff overestimates the ef- 
fective charge of the polyelectrolyte 
counterion complex. 

The effective friction of the poly- 
electrolyte counterion complex is in 
agreement with the hydrodynamic 
friction (r = 1/D) for short chains 
only. For longer chains beyond the 
Debye length in the system (indi- 
cated by the dashed lines) the effec- 
tive friction deviates and seems to 
approach a linear behaviour. 

The effective charge Q e s and the 
effective friction T e fj per monomer 
show a strong dependence on the 
chain length N for short chains. For 
long chains, both values are con- 
stant. The transition occurs on 
a length-scale similar to the Debye 
length (dashed lines). 

The interplay between counterion 
condensation and conformational 
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change of the polyelectrolyte results 
in the combined behaviour of effec- 
tive friction and effective charge that 
leads to the observed length- depen- 
dence of the electrophoretic mobility. 
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Figure 1: Both the end-to-end distance R e and 
the radius of gyration R g at vanishing exter- 
nal field for a monomer concentration of Cpe = 
Cci = 10 niM exhibit an effective scaling ex- 
ponent v = 0.85. The hydrodynamic radius 
i?h shows a different behaviour for short chains 
and only slowly reaches the asymptotic scaling. 
The static chain properties are not influenced by 
hydrodynamic interactions (HI) and fully agree 
with the ones obtained without HI. 
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Figure 2: (a) The end-to-end distance R e and 
(b) the radius of gyration radius R g for a chain 
of N = 40 monomers differ from the value at 
zero external field E (dashed line), if a criti- 
cal value is reached. Beyond this threshold, the 
polyelectrolyte chain is extended (increased R e 
and i?g) and starts to align with the field di- 
rected in x-direction. This effect is independent 
of hydrodynamic interactions (HI) . 
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Figure 3: The end-to-end distance R e and the 
radius of gyration R g at zero external field and 
for a monomer concentration Cpe = 10 mM and 
no added salt (c s = M) exhibit an effective 
scaling exponent v = 0.85 (dashed line). With 
c s = 1 M monovalent salt added, the scaling 
exponent drops down to v — 0.68 (dotted line). 



23 



1 



0.9 




0.8 



> 



0.7 



0.6 



0.5 



0.1 



-l 



Figure 4: The effective scaling factor v shows 
strong dependence of the Debye length. With 
out added salt (i.e. large Debye lengths) v is 
close to 1, whereas for high salt concentrations 
(i.e. small Debye lengths) approaches the Flory 
number, the value of an uncharged polymer. 
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Figure 5: As a simple estimate of the number 
of condensed counterions iVci all ions within 
a distance of 2cr around the PE chain are 
summed up. For long polyelectrolyte chains, the 
value predicted by Equation [T] is approached. 
The static counterion distributions are not in- 
fluenced by hydrodynamic interactions (HI). 
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Figure 6: The number of condensed counteri- 
ons iVci for a polyelectrolyte chain of = 40 
depends on the applied external electrical field. 
Above a threshold of E — 0.1, counterions start 
to be stripped from the PE chain and Nqi dif- 
fers from the value at zero field (dashed line). 
This effect is independent of hydrodynamic in- 
teractions (HI). 
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Figure 7: The diffusion coefficient D of a poly- 
electrolyte chain of length N = 32 is deter- 
mined, (a) A fit (dashed line) to the linear part 
of the mean square displacement yields a dif- 
fusion coefficient of D = 0.0045 ± 0.0002 via 
Equation [2j (b) By fitting a t~ 3 ^ 2 power law to 
the long-time tail of the center of mass auto- 
correlation function and using Equation [3] we 
obtain D = 0.0041 ± 0.0005. 
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Figure 8: The diffusion of a polyelectrolyte 
chain of length N at a monomer concentration 
c m = lOmM, normalised by the monomer diffu- 
sion D Q , shows the influence of hydrodynamic 
interactions (HI). With HI a scaling exponent 
of m = 0.63 ± 0.01 is obtained (solid line), 
whereas without HI, the slope is m = 1.02±0.02 
(dashed line). The diffusion in presence of HI 
agrees with the values obtained from Equa- 
tion HI where Rh is determined from the sim- 
ulation (triangles). 
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Figure 9: The translational diffusion coefficient 
measured at a reduced electric field E = 0.1 
(circles) is in agreement with the measurements 
at zero electric field (solid line). 
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Figure 10: The normalised electrophoretic mo- 
bility at a monomer concentration of c m = 5 
mM with hydro dynamic interactions (HI) repro- 
duces the experimentally observed behaviour 
and shows a maximum for intermediate chains. 
Without HI the measures mobility strongly de- 
viates from this showing a decrease with increas- 
ing chain length N. This behaviour can be ex- 
plained by Equation [7| in the local force picture 
(dashed line) and approaches the limiting value 
of the Manning prediction (Eq. [5]). 
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Figure 11: Electrophoretic mobility /x of poly- 
electrolyte chains at monomer concentration 
c m = 5 mM measured without electric field (cir- 
cles) and with field at E = 0.05 and 0.1 (trian- 
gles) with the corresponding error bar. (Data 
sets have been shifted to increase readability.) 
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Figure 12: The normalised electrophoretic mo- 
bility h/hq shows a strong dependence on the 
monomer concentration c m . For dilute systems, 
the experimental behaviour is recovered, while 
at high monomer concentrations, the shape is 
significantly altered. 
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Figure 13: The approximate position of the 
maximum N max changes with the Debye length 
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Figure 14: The integrated ion distribution 1(d) 
for a chain of iV = 32 monomers shows a clear 
inflection point identified by the zero of the sec- 
ond derivative with respect to hx(d). The value 
of 1(d) at the inflection point specifies the frac- 
tion of condensed counterions fci- 
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Figure 15: Estimated charge Q^g of polyelec- 
trolyte chains of varying length N using Equa- 
tion [12] with and without hydrodynamic inter- 
actions (HI). 
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Figure 16: The ion velocity distribution fci(^) 
for a polyelectrolyte chain of length N = 32 
shows the transition from co-moving counteri- 
ons with a velocity similar to vpe (dashed line 
and dotted line indicate vpe with and without 
hydrodynamic interactions (HI)) to the free ion 
velocity vqi,o (dash-dot line). vci(d = do) = 
defines the threshold between co- and counter 
moving ions and is used to determine the ef- 
fective charge of the polyelectrolyte counterion 
complex in Equation [TU 
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Figure 17: The effective charge of polyelec- 
trolyte chains of length N using the estima- 
tors Qeff t° Qeff- The dotted line indicates 
the bare, unscreened charge of the polyelec- 
trolyte, whereas the solid line shows a predic- 
tion based on counterion condensation theory, 
with £ = 2.84 being the condensation parame- 
ter. Whereas estimators to agree over 
the range of lengths iV tested, overesti- 
mates the effective charge of the polyelectrolyte 
counterion complex. 
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Figure 18: The effective friction of the polyelec- 
trolyte counterion complex is in agreement with 
the hydrodynamic friction (T — 1/D) for short 
chains only. For longer chains beyond the Debye 
length in the system (indicated by the dashed 
lines) the effective friction deviates and seems 
to approach a linear behaviour. 
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Figure 19: The effective charge Q e g and the ef- 
fective friction F e Q per monomer show a strong 
dependence on the chain length iV for short 
chains. For long chains, both values are con- 
stant. The transition occurs on a length-scale 
similar to the Debye length (dashed lines). 
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Figure 20: The interplay between counterion 
condensation and conformational change of the 
polyelectrolyte results in the combined be- 
haviour of effective friction and effective charge 
that leads to the observed length- dependence 
of the electrophoretic mobility 
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